Neutrinos as a Diagnostic of High Energy Astrophysical Processes 
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A leading candidate for the extragalactic source of high energy cosmic rays is the Fermi engine 
mechanism, in which protons confined by magnetic fields are accelerated to very high energy through 
repeated scattering by plasma shock fronts. In the process of acceleration, collisions of trapped pro- 
tons with the ambient plasma produce pions which decay to electromagnetic energy and neutrinos. 
For optically thin sources, a strong connection between the emerging cosmic rays and secondary neu- 
trinos can be established. In this context, we show the feasibility of using the Glashow resonance 
as a discriminator between the pp and py interactions in Fermi engines as sources of neutrinos. In 
particular, we demonstrate how three years of observation at the km 3 IceCube facility can serve as 
a filter for the dominance of the pp interaction at the source. 



Neutrinos can serve as unique astronomical messen- 
gers. Except for oscillations induced by transit in a vac- 
uum higgs field, neutrinos propagate without interactions 
between source and Earth, providing powerful probes of 
high energy astrophysics [lj • The deployment of the km 3 
IceCube facility at the South Pole @ will greatly increase 
the statistics required for the realization of such a pro- 
gram. 

The flavor composition of neutrinos originating at as- 
trophysical sources can also serve as a probe of new 
physics in the electroweak sector Specifically, Ice- 
Cube will have the capability to clearly identify neutrino 
species [U, and consequently will be able to measure 
deviations of flavor composition from standard expecta- 
tions. In addition, in resonant scattering v e e~ — > W~ — > 
anything p|, the detector can simultaneously discrimi- 
nate between is e and V e . In this Letter, we show that the 
signal for V e at the Glashow resonance, when normalized 
to the total v+v flux, can be used to differentiate between 
the two primary candidates (jr/ and pp collisions) for 
neutrino-producing interactions in optically thin sources 
of cosmic rays. 

It is helpful to envision the cosmic ray engines as ma- 
chines where protons are accelerated and (possibly) per- 
manently confined by the magnetic fields of the accelera- 
tion region. The production of neutrons and charged pi- 
ons and subsequent decay produces both neutrinos and 



cosmic rays: the former via 7r + 



e^v e v,,v. 



(and the 



conjugate process - more on this below), the latter via 
neutron diffusion from the region of the confined pro- 
tons. If the neutrino-emitting source also produces high 
energy cosmic rays, then pion production must be the 
principal agent for the high energy cutoff on the pro- 
ton spectrum 0, Q • Conversely, since the protons must 
undergo sufficient acceleration, inelastic pion production 
needs to be small below the cutoff energy; consequently, 
the plasma must be optically thin. Since the interac- 
tion time for protons is greatly increased over that of 
neutrons because of magnetic confinement, the neutrons 
escape before interacting, and on decay give rise to the 
observed cosmic ray flux. A desirable property of this 
low-damping scenario is that a single source will produce 
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cosmic rays with a smooth spectrum across a wide range 
of energy. 

For optically thin sources, the neutrino power densit 
scales linearly with the cosmic ray power density e CH 
The actual value of the neutrino flux depends on what 
fraction of the proton energy is converted to pions (which 
then decay to neutrinos and photons). To quantify this, 
here we define as the ratio of pion energy to the emerg- 
ing nucleon energy at the source. Since there is about one 
cosmic ray (neutron) produced per proton collision, a sig- 
nificant consequence of the neutron leakage model is that 

is expected to be approximately equal to the fraction 
of collision energy carried off by pions in a single colli- 
sion. Following Waxman and Bahcall [g, the expected 
total (all species) neutrino flux is 

(1) 

where for no source evolution £ 2 w 0.6, whereas for an 
evolution oc (1 + z) 3 as seen in the star-formation rate, 

Depending on the relative ambient gas and photon 
densities, pion production proceeds either through in- 
elastic pp scattering |Tcj . or photopion production pre- 
dominantly through the resonant process P7 — » A + — > 
nTT + OTpTr° For the first of these, pion production 
is well-characterized by a central rapidity (Feynman) 
plateau, yielding sa 0.6 11]. For resonant photo- 
production, the inelasticity is kinematically determined 
by requiring equal boosts for the decay products of the 
A+ [H, giving e T « 0.25 0. _ 

The first production mechanism (described above) in- 
volves the dominance of inelastic pp collisions in gen- 
crating charged pions. The nearly isotopically neutral 
mix of pions will create on decay a neutrino population 
in the ratio N v — Np = 2N Ve = 2Np e . In propa- 
gation to Earth, a distance longer than all oscillation 
lengths, flavor-changing amplitudes are replaced by prob- 
abilities , resulting in equal fluxes for all six neutrino 
species, so that jvi? arth = | N^~. Here, we have incor- 
porated maximal <-» v T mixing and the known small- 
ness of |(zv e |f3)| 2 , where ~ (v^ + z/ T )/V2 is the third 
neutrino eigenstate. 
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The second mechanism proceeds via photoproduction 
of pions by trapped protons on the thermal photon 
background, leaving the isotopically asymmetric process 
P7 — * A + — * 7r + n as the dominant source of neutrinos. 
At production, N Vli = N Vpi = N Vc > N Ve . After oscilla- 
tion, the V e flux at Earth is 



N-, 



Earth 



- sm 2 d Q COS 2 6 G A^total 



(2) 



Using the most recent SNO value of the solar mixing 
angle Q ~ 32.5° (H, we obtain N^ arth = i N^. 
Because of averaging over various thermal environments, 
this relation holds across the neutrino spectrum || . 

The preceding discussion indicates the feasibility of 
differentiating between the two processes as follows: (i) 
obtain a normalization for the total flux from data out- 
side the Glashow resonance region; (ii) measure the event 
rate in the resonance region, and compare with the ratios 
given above. In what follows, we will assess the potential 
of IceCube to obtain a statistically significant signal-to- 
background ratio. 

The cross section for Glashow resonant scattering reads 



there will be a single electromagnetic burst triggered by 
the electron interaction or by the r decay, each however 
manifesting only half the energy. As a conservative proce- 
dure, we will use in our signal estimate only the hadronic 
events with full energy in the resonant band, giving (in 
the case of pp interactions) a signal dAf jdt — 3.2 yr _1 . 
The background will then constitute the non-resonant 
scattering of v e and V e within the resonant band. This is 
calculated as an integral of the flux times the charged cur- 
rent cross section [l7j over the resonant acceptance bin 
(10 6 - 7 GeV, 10 6 ' 9 GcV) of the AMANDA detector 
The process being off-resonance, the Earth is transparent 
and so we take Ail = An. The number of target nucleons 
is twice the number of electrons. This procedure yields 
a background of 0.6 events yr -1 . Thus, within 3 years 
of data accumulation (1/5 of the total lifetime of the ex- 
periment), we expect about 2 background events, with a 
signal of 9.6 events, well above the 99%CL for discovery 
(corresponding to 6.69 events (l9jl- 

Electron antineutrinos can also be produced through 
neutron /3-decay. This contribution turns out to be neg- 
ligible. To obtain an estimate, we sum over the neutron- 
emitting sources out to the edge of the universe at a 
distance 1/H |2C| 



|9Jl| 2 5{2m e E v - m 2 w ) 



(3) 



2m e E v 

where \Wl\ 2 = g 2 m 2 v /2, m e and mw are the electron and 
W masses, and g 2 — 4ma{mw)/ sin 2 0w(jnw) — 0.43. 
Now, for an incoming flux of antineutrinos J-jy c (E v ) the 
event rate (assuming 100% detection efficiency) reads 

dAf 



— = N cS Ail J dE v J Vs {E v ) a{E u ) 
4m e V 2m e 



(4) 



where N c g ~ 6 x 10 38 is the number of target electrons 
in the effective volume V ~ 2 km 3 of the IceCube ex- 
periment, and Ail ~ 2ir is the solid angle aperture |16| . 
If neutrinos are produced via pp inelastic collisions, ac- 
cording to the previous discussion (with = 0.6 and 
£ z — 3), the antineutrino flux on Earth is E 2 Jv e {E v ) = 
6 x 10~ 9 GeV cm~ 2 s _1 sr _1 , so that the event rate is 
dAf jdt = 4.6 yr . For antineutrinos produced in pho- 
topion interactions (e^ = 0.25), the expected rate is a 
factor of 6 smaller, dAf jdt = 0.8 yr" 1 . 

Can the expected resonant signal be differentiated 
from the continuum background? In order to estimate 
the background, we note that all but the muon decays of 
the W devolve into an electromagnetic shower in the de- 
tector. For the hadronic decays (which constitute 70% of 
the total), the interaction mean free path of charged pions 
in ice is orders of magnitude smaller than the pion decay 
length at TeV energies, allowing a progressive channel- 
ing of nearly all the energy into electromagnetic modes 
through 7r° decay. For the e~v e and t~v t decay modes, 
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(5) 



where T n (E n ) is the neutron volume cmissivity and m n 
the neutron mass. Here, we have assumed that the neu- 
trino is produced monoenergetically in the neutron rest 
frame, i.e., e ~ Smjf{l — m 2 /8 2 mM)/2 ~ 0.55 MeV, 
where 5m n ~ 1.30 MeV is the neutron proton mass dif- 
ference. An upper limit can be placed on T n by assum- 
ing an extreme case in which all the cosmic rays escaping 
the source are neutrons, i.e., e CR = J dE n E n T n (E n ). 
With the production rate of ultrahigh energy protons 



;[10 10 ,10 12 ] 



5 x 10 44 erg Mpc -3 yr" 1 |2l|, and an as- 



sumed injection spectrum T n oc E n , we obtain 
E 2 $ F « 3 x 10~ n GeVcm^s^sr- 1 , 



(6) 



about 3 orders of magnitude smaller than the charge pion 
contribution with = 0.6. (Note that oscillations will 
reduce the v e flux on Earth by 40% from this value.) 

Finally, we note that in the case that pions are 
produced via pry — > A + — > V7r, the event rate for 
the Glashow resonance cannot be separated from back- 
ground. Thus, the presence or absence of a signal at 
the Glashow resonance can be used as a filter for pp- 
dominance at the source. 

In summary, we have shown how several years of 
data collection at IceCube can potentially isolate the 
dominance of inelastic pp interactions as neutrino 
progenitors in cosmic ray sources. The analysis depends 
on the sources being optically thin, which can be 
ascertained if the observed diffuse neutrino flux is in 
agreement with Eq. (JTJ- Then a signal at the Glashow 
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resonance will provide a strong marker for pp-dominance. 

Note added: The previous version of this paper 
contained a brief discussion of how the flux of neutrinos 
from optically thin sources can be used as a marker 
of the cosmic ray Galactic/extra- galactic transition. 
This topic has received detailed discussion in a separate 
paper |22j |. 
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